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Abstract 

In the Principia Mathematica Sir Isaac Newton proved that concentric mass shells with equivalent 
mass distributions produce the same external gravitational field and thus that the problem of estimating a 
continuous interior mass distribution from external field information alone is ill-posed. What is generally 
less well known is that finite collections of point masses contained in some bounded domain produce 
a unique field in the exterior domain, which means that the associated basis functions (often called 
"fundamental solutions") are independent. A new proof of this result is given in this paper that can be 
generalized to other finite combinations of point source distributions. For example, one result this paper 
shows in is that a finite combination of (gravitational or electrostatic) point dipole sources contained 
in some interior region produces a unique field in the corresponding exterior region of interest. 

Since no direct proofs of uniqueness results of this type are known for the setting, which is the 
setting of primary practical interest, an indirect strategy is necessary. The strategy employed in the 
paper is to develop results for analytic functions in the complex plane, C, where logarithmic source basis 
functions correspond to point mass basis functions, and then carry them over to harmonic functions in 
the real plane, R^, and from there to harmonic functions in M^. Although some of the results obtained 
can be generalized from to M", for n > 3, far more results are shown for and for C than are 
shown for these more general settings. For example, in the complex plane, the paper shows that a finite 
combination of higher order poles of any order in the interior of a unit disk always corresponds to a 
unique analytic function in the exterior of a unit disk. 

Key words: Laplace's equation, inverse problem, potential theory, point sources, Vandermode 
matrix, fundamental solutions, multipole 
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1 Applied Backdrop 

Inverse source problems associated with harmonic functions (i.e., ones that satisfy Laplace's equation) are 
a small area of modern inverse source theory; however, this area encompasses problems in geophysics, 
geoexploration and electrostatics, as well as many other applied sciences. Moreover this area undoubtedly 
contains the oldest substantial mathematical result of any real significance associated with inverse source 
theory: Sir Isaac Newton showed in his Principia Mathematica that a continuous spherically symmetric 
distribution produces an external field equivalent to that of a point mass located at the sphere's center, 
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provided that the total masses of both are the same. From this resuh, it follows immediately that the 
problem of determining continuous mass distributions inside some interior region from external gravity field 
information alone is an ill-posed problem, because many solutions are possible. What is generally less well 
known is a proof exists for R" (n > 2) that shows that a finite set of point masses located inside some 
bounded region produces a unique field in the exterior region [f 2 . This paper presents an alternative proof 
of this R" result, as well as various generalizations of it. One of these generalizations is that a finite set of 
point dipoles is shown to be independent in R'^. Other generalizations of this result are shown in the complex 
plane, C, and the real plane, R^. Finally, one open question is analyzed extensively here; namely, are linear 
combinations of R'^ combined point mass/point dipole basis functions linearly independent? In the end, the 
independence of these basis functions is shown to hinge on the invertibility of a matrix with complex entries 
that is related to the Vandermode matrix [and which is specified by (|59p]. Clearly, in general, this generalized 
Vandermode matrix is invertible, however, currently no proof of this invertibility for all n is known. 

In what follows, the term point source will be used to denote a general source term [i.e., a point mass 
or (point charge), a point mass dipole (or an electrostatic dipole), a point quadrapole or a higher order 
(electrostatic) multipole in M" (n > 2) — or a logarithmic pole (with a specific form), a simple pole or a 
higher order pole in C]. 

Because uniqueness results may seem to be somewhat removed from the applications arena, before con- 
sidering the paper's basic approach and content in more detail, it is appropriate to briefly consider the 
motivations for studying point source uniqueness results. The primary motivations can be framed as follows: 

1. By better understanding and quantifying various point source uniqueness results, the hope is that 
geophysical inverse source theory can be better understood and placed on a firmer physical and math- 
ematical foundation. (Point source uniqueness results should be viewed as merely a first tentative step 
in this goal.) 

2. Point source uniqueness results show that, in theory, when well formulated algorithms are properly 
implemented, point source determination software always produces reliable results. In particular, 
by showing that combinations of point sources and point dipoles are linearly independent, one can 
immediately infer that the DIDACKS implementation that interpolate for the scalar potential of gravity 
and the vector components of gravity are mathematically well framed, as discussed in [6l [7l [8] and [9] . 

3. As discussed in Appendix A of [8j, this uniqueness result for point sources and point dipoles also 
shows that geophysical collocation procedures for similar data sets (such as geoid height and the 
vector components of gravity disturbance) are mathematically consistent since the associated error- 
free covariance matrix can always be inverted. 

4. Finally, as always, when a different topic is addressed, different mathematical techniques may be called 
for, and thus there is the possibility of uncovering new mathematical theorems, results and techniques 
that may be of interest in other domains. Here, not only are the results dealing with logs and poles in C 
suggestive, but so are several of the (geometric) theorems that allow for results to be carried over from 
R^ to R". Also, since Vandermode's matrix in an extended form is used here, it is not unreasonable to 
think that there are mathematical links to general interpolation theory, where the Vandermode matrix 
is often employed, and that this is a two-way street (which means that there may be general theorems 
in interpolation theory that allow for more general point source uniqueness theorems than those proven 
here) . 

Here item I. is worth delving into a little deeper, since placing inverse source gravity theory (i.e., R'^ Laplacian 
inverse source theory) on a firmer mathematical and physical foundation is one goal of the DIDACKS 
sequence of papers initiated by [7]. In this regard, there are a number of interconnected points that are 
relevant, and it is appropriate to point one of them out here. First, consider the point made in 4 that since 



2 



each point source can be replaced by an equivalent sphere of uniform density, one can replace a distribution 
of point sources with an equivalent set of spheres that produces a more uniform distribution of mass. Thus, 
consider a collection of point masses that have different strengths from point to point and are located on a 
regular cubic grid. Since the grid spacing is uniform each point mass can be replaced by a sphere of like size 
that has a corresponding density. From the uniqueness theorems given here one can infer that the estimation 
problem for the associated densities of this spherical arrangement is well posed. Moreover, since it is well 
known that as the grid spacing for a set of point masses increases the underlying estimation problem becomes 
better conditioned, one can immediately infer that the same holds for this spherical approximation. This 
result has other generalizations and interpretations as well. For example, in geoexploration it is common 
to use a set of right parallelepipeds of uniform density to approximate continuous mass densities since they 
produce potentials that can be calculated in closed form [51 p. 398] and they yield inverse source density 
estimation equations that are (relatively) stable. As a next step, one might consider extending the uniqueness 
results here to uniform cubic latices (or one might consider applying DIDACKS theory to the estimation of 
such cubic distributions) . This and other points raised above clearly dovetail with similar discussion in '10] . 



2 Mathematical Preliminaries 

The general method of proof employed is to start with the easiest to handle uniqueness case — which turns 
out to be a distribution of simple poles of the form l/{z — Zk) in the complex plane-and then generalize this 
result in various ways: First to uniqueness of logarithmic basis functions in C, then to higher order poles of 
the form l/(z — Zfc)". Next a (well-known) one-to-one correspondence between various types of expansions 
in C and is noted: logarithmic potentials point masses; simple poles dipoles; poles of order 

n multipoles of higher order for n = 1, 2, 3. It is then shown that uniqueness of an expansion in C 

implies uniqueness of a like expansion in — thus proving the cases. Finally it is shown explicitly that 
point mass uniqueness results in R^ imply point mass uniqueness results in M'^. These correspondences and 
the symbols used to represent the associated expansions in the sequel are shown in Table [TJ The arguments 
employed can clearly be generalized in two ways: (1) To show that similar results hold for multipoles of 
order three or less in R^. (2) To show that results obtained in imply similar results in R^ for N > 2. 
For concreteness and for ease of exposition these two generalizations will not be explicitly dealt with here 
since they do not correspond to commonly occuring physically relevant cases and, when handled explicitly, 
entail some technical difficulties. 

In R^ and R"^ the conservative vector fields of interest can be obtained by taking the gradient of a scalar 
field: 

FiX) = -VUiX), (1) 

where, of course, V is the iV— dimensional gradient and where X G R^. Potentials specified by ^ satisfy 
Laplace's equation in TV— dimensions: V^J7 — 0. Arguably the most significant cases, and the ones that will 
be focused on here, are the electrostatic and gravitational fields in R^ and R"^. Because of the issue of the 
sign associated with the mutual attraction of point sources and the issue of the sign of the energy density 
play no role in what follows, the only relevant difference between the gravitational and electrostatic cases 
is that point mass strengths (nik) are generally assumed to be non-negative, but electric change strengths 
(qk) can have either sign. Since point source uniqueness results that hold for sources of either sign also 
clearly hold when the sources are all assumed to be positive and since a historical precedent exists for stating 
uniqueness results in terms of point masses, in the sequel it will always be assumed that while the scalar 
point source parameters are called a point masses and denoted m^, their assigned values can take on either 
sign. Likewise the vector point source distribution that corresponds to a electrostatic dipole or point mass 
dipole will be denoted Dk and it will be called a point dipole. Higher order point multipoles will also be 
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Description 


Symbol 


C Analog 


C Symbol 


Harmonic Functions 


U 


Analytic Functions 


/ 


Point Masses 


V 


Logarithm Terms 


9 


Point Dipoles 


W 


Simple Poles 


h 


Multipoles 




Higher Order Poles 





Table 1: Corresponding Types of Potential and Analytic Functions 



considered in what follows and point masses, point dipoles and point multipoles will be collectively referred 
to as point sources. In all cases, since either sign is allowed for these point sources and since uniqueness 
results pertain to whether the source terms form a linearly dependent set of basis functions in the exterior 
region of interest, the overall choice of sign convention for these basis functions and associated point source 
strengths does not matter here and the basis sign conventions are chosen for convenience. As previously 
noted, while point source distributions for for iV > 3 exist and the results here can be naturally extended 
to handle them they will not be explicitly considered here. Finally, observe that magnetic dipoles have a 
different form from the gravitational or electrostatic (and magnetostatic) potential forms assumed here and 
the question of adapting the results here to their study will also not be addressed. 

First, as a general convention, let Nk denote a finite integer greater than zero and let the subscript k (or 
k') be used to index the Nk point sources under consideration so that k and k' — 1, 2, 3, . . . , Nk is always 
understood. To further fix notation in for > 2, let the Nk fixed distinct point-source locations be 
specified by Xk so that Xk Xk' when k' ^ k. It will also be assumed that all of the sources are located in 
some bounded domain. 

Introducing a specific symbol for point mass potentials, V , instead of the general symbol U used in ([1]) 
results in the following definition 

Nk 

ViX) = J2 "^fc In i\X - Xk\-') (2) 



fe=i 



m 



m 



and 



iVj 



R^ 

Likewise the point dipole potential form is 



W{X) = Y, Dk-V lii{\X - Xk\-') (4) 



k=l 

in R^ and 

Nk 

WiX) = Dk-V i\X ~ Xk\-') (5) 

k=l 

in R'^. The relationships of these potentials and their corresponding complex analytic counter parts are 
shown in Table [1] These correspondences will be addressed below. 
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Next consider uniqueness results stated in terms of scalar potentials. So long as Laplace's equation is 
satisfied in the region of interest the exact shape of the exterior region is immaterial due to the uniqueness 
of Dirichlet boundary value problems so, without loss of generality, assume that the mass distributions are 
located in some bounding sphere. Further, since the origin of coordinates and the length scale also does 
not change, for the desired final uniqueness results it can be assumed that the harmonic region of interest 
is \X\ > 1 and that the sources are in the compliment of this region: < \Xk\ < 1 (where, for later 
convenience, it is assumed that the origin is not situated directly over any particular source). Point mass 
uniqueness requirements can then be stated as the condition that if V{X) = 0, for all X > 1, then nik = 
for all k and that the converse also holds. Which is to say that for any finite Nk > 0, if rrik ^ 0, for all fc, 
then V(X) ^ for some values oi X > 1. Dipole uniqueness can be stated similarly: If V{X) — 0, for all 
\X\ > 1, then Dk = for all k, or conversely as the condition that for any finite N/, > if D^. ^ 0, for all 
fc, then V(X) ^ 0. It is assumed throughout that any nonzero values of ruk and Dk are bounded (while the 
fixed nature of X^ rules out the formation of dipoles in the limit rn^ — > oo and mk' — > — oo with \mk\ = \mk> \ 
for some pair of points indexed by k and fc', the unbounded mass case is still best bypassed). 

Several observations are relevant. By introducing the basis functions 

■^k{X) - In ^ ^ ^ - In ^ = In ^*^L , (6) 
\X-Xk\ \X\ \X-Xk\ 

Equation ^ can be reexpressed as 

V{X)^Y. rnk^kiX) +moln ^ 
fc=i ^ 

where rrio = X^aSi '^fc- The case toq can be easily disposed of since — > as \X\ — > oo and it is clear 
that |y| — > oo as |X| oo unless mo = 0. Consequently, without loss of generality take mo = 0, so that the 
form 

F(X) mfc*fc(X) (7) 
fc=i 

is always assumed for in what follows. Likewise all of the harmonic functions considered here will be 
assumed to vanish at infinity by convention. 

In R^ , uniqueness results for a scalar potential imply uniqueness results for the vector field itself as can 
be seen from the following argument. Consider the line integral 

U{X)^U{Xo)+ I yU-di (8) 

where I — £{s) denotes a parameterized path (with arclength s). The contention is that since X and Xq 
are arbritrary points in the exterior region and £(s) is also assumed to lie wholly in this exterior region, 
F ^ ■^=^ V ^ and conversely F — <=4> U — 0, where F and U are related by ([1]). For example, 
if F{X') < at some point X' > 1, then due to the mean value theorem for harmonic functions VC/ > 
holds in some finite neighborhood of X' so that both X and Xq, along with the path connecting them in ([8]), 
can be taken to be inside this same neighborhood. This, in turn, means that at least one of the two values 
U{X) or U{Xo) must be nonzero. Alternatively, if U{X) > 0, then Xo can be set to a point at infinity and 
from ([S]) it is clear that F ^ must occur at someplace along the line integral. Without loss of generality, 
uniqueness results will thus be stated in terms of scalar potentials for convenience. 
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Uniqueness in the complex setting is addressed first since this is the easiest route to the desired results, 
which can be readily generalized to K^. It is useful to have a common (and commonly used) symbolism for 
addressing uniqueness issues in and C. Let x and y denote standard Cartesian coordinates in either 
setting: In M^, X = {x, y)^ and = {xk, yk)"^ (where T denotes the transpose); while in C, z = x + iy 
and Zk — Xk + i yk- In both settings y^x^ + y^ > 1 and y/ xj, + < 1. Recall from elementary treatments of 
analytic functions that there is a general mapping between harmonic functions defined over some subregion 
of and analytic functions defined over the corresponding subregion of the complex plane as indicated by 
Table [TJ This mapping can be done uniquely when certain reasonable conditions are met with regards to 
branch-cuts and the nature of the region under consideration and it is assumed that the reader is familiar 
with them. Specifically, if U is harmonic in let f{z) denote the unique analytic function in C whose real 
component corresponds to C/(x, y), in which case f{z) will be called the standard completion of U{x, y). 
The standard completion of the sum of real functions obeys the principle of linear superposition, so that, for 
example, the standard completion of is a linear superposition of terms that are the standard completions 
of the ^'fc(X). Specifically, let standard completions of ^'^(X) be denoted ipk{z)- Since the real part of In 2; 
is In |z|, it is obvious from ((6|) that the standard completion of '^k is given by 

V'fe(z) = In . ^ . In — \n - (9) 

(Z- Zkj Z- Zk z 

and thus (using the notation indicated in Table [1] for the corresponding complex logarithmic case) g{z) can 
be written as 

.9(2) = ^Mfe V'fe(^) , (10) 

fe=l 

where /ik & C 

As before, a general series based on In [1/ z — Zk] can be contemplated here since it corresponds to adding 
a term involving 1/z to (jlOp : however, it is readily proved that when this term is present it always causes 
g{z) ^ for large z and so its presence need not be considered further. This might all seem straightforward, 
but even here some care is called for. Moreover, since the basic strategy used in the sequel will be to prove 
uniqueness (i.e., linear independence) in one setting and then to obtain uniqueness results in all other settings 
by applying a linear uniqueness preserving mapping, the general theorems that are required are best stated 
explicity (for clarity and uniformity of exposition). Even though the underlying concepts are well known, 
these linear independence preserving theorems do not necessarily take conventional forms. 



3 Uniqueness Preserving Mappings 

This section discusses uniqueness preserving mappings that are used in the sequel. As noted at the end of 
the last section, the first of these mappings is associated with the act of standard analytic completion and 
entails a unique correspondence between analytic functions in the complex pane and harmonic functions in 
M^. Since uniqueness results will first be shown in the complex setting and then mapped into the setting, 
first consider what uniqueness means in the complex setting: 

Definition [3jl A set of basis functions {fk}kli is said to produce a unique expansion in the complex 
setting when they meet the following criteria. First each basis function must be a bounded analytic function 
for 1 2; I > 1. Second each basis function must vanish at infinity. Third, a sum of the form 

f{z) = '^fikfk{z) with fik & C (11) 

fc=i 
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must be linearly independent for \z\ > 1, where linear independence means that f{z) = holds if and only 
if /I A; = for all k. 

Let Re{/} denote the real part of / and Im{/} the imaginary part, then if Uk{x, y) = Re{fk{z)}, 
Vkix, y) = Im{/fe(z)}, ak = Re{fik} and (3k = lm{fik}; from fk{z) = Uk{x, y) + ivk{x, y) and ^fe = ak+iPk 
it follows that 

^(^{f} ^^^(^{t^klkiz)} ^^[akUk{x, y) - (3kVk{x, y)]. (12) 

k=l k=l 

Since /ifc = for all k implies that both = and = hold for all k, uniqueness of the analytic set 
of basis functions {fk}i '' , implies simultaneous uniqueness of the pair of conjugate harmonic basis function 
sets {uk{x, y)}k=i and {vk{x, y)}k=i- Uniqueness in the setting is defined analogously to Definition [Sjl. 
This result can be summarized as a theorem: 

Theorem 1. // the set of basis functions {fk}k=i '^'^^ unique in the complex setting and if Uk{x, y) = 

; tl 

are linearly independent or unique in the setting. 



Re{fk{z)} andvk{x, y) = Im{fk{z)} , then the combined set of basis functions {uk{x, y)}^Jli^{vk{x, y)}^^i 



Notice that since the act of standard harmonic completion is unique, if uniqueness can be shown in for 
a sequence of harmonic conjugate pairs then uniqueness for basis functions of the form (jlip follows. 

For the complex setting, in what follows only two general types of basis functions will need to be consid- 
ered: (1) Logarithmic potentials discussed above of the form ipk{z) (which, as discussed latter, are analytic 
from the branch-cut considerations in the first part of Appendix A). (2) Poles of the form jik/iz — Zk)", for 
finite n > 0. For f{z) = 'Yl!k=il^k4'k{z), the restriction Pk = can be made since the argument dependent 
parts of the logarithmic term occuring in ([9]) [i.e., Im {■(/'^(z)}] are not of general interest in the harmonic 
setting. Uniqueness results in the complex plane for a series of logarithmic basis functions can thus be used 
to show point mass uniqueness in M^. 

Next consider a series of poles of fixed order n: 

k=i 

The case n = 1, h{z) = h^^^{z), is of special interest and results in the well-known linear combination of 
simple poles: 

Uniqueness results will first be shown for an expansion of this form. Since 

Mfc _ fikiz* - zl) _ ak{x - Xk) + (3k{y - yk) , . Pkjx - Xk) - akjy " Vk) 



z-zk \z-zk\'^ {x - XkY + {y - ykY {x - XkY + {y ~ ykf 

the real part of h{z) corresponds to an dipole expansion W{X) given by ([l]), as will now be explicitly 
shown. 

The various point dipole terms contained in Q are also known as first order multipoles. For future 
reference it is useful to have a consistent notation for delineating multipole basis functions of various orders. 
The order of a multipole corresponds to the number of subscripts it has, so a first order multipole has a 
single subscript that can take on the values 1, 2, 3, . . . , in M^. Thus the subscripts of a dipole (or first 



7 



order multipole) basis function are associated with the various directions in R . In R this subscript takes 



on two values so that the two first order multipole basis functions associated with the source position 
Xk can be written as M^''' = M.j{X, Xk) for j = 1, 2. The R^ dipole or first order multipole basis function 
oriented along the a;-axis in is given by 

- ^ In ^===== = -z:— T^^ri— ■ (16) 



I 








X - 


- Xk 




y/{x - Xfe)2 - 




- VkY 


[x- 


xkY 


+ {y 


-VkY 


aligned along 


the 


y-axis in 


is I 


^iven 


by 




1 










- Vk 




^{x - XkY H 




- VkY 


(x^ 


xkY 


+ {v 


^VkY 



Thus uniqueness results in the complex plane for a series of simple poles can be used to show dipole uniqueness 
in the R^ setting. Theorem [1] (and its reverse) clearly involves a change of setting form C to R^ (or R^ to C) 



4 Uniqueness of Complex Poles of The Same Order 

This section considers relatively straightforward results about uniqueness of expansions of poles and unique- 
ness of expansions of logarithmic basis functions. In particular, the theorem that shows that an expansion 
in terms of simple poles, h{z) as given by (I14p . is unique is readily stated and proved. The analogous result 
is also easy to prove for expansions of higher-order poles of a given type (i.e., where all the poles are of some 
specified order), as well logarithmic basis functions. A consideration of the more difficult uniqueness results 
for mixed types of expansions are postponed until Section [7] and, in the end, concrete proofs of these mixed 
uniqueness results prove to be elusive. 

Before proceeding to a statement of the desired theorem, several observations are in order. First, if 
2fc = for any k, then it is simply necessary to translate and rescale, when required, so that ^ can 
be assumed. Second, as in R"^ for concentric spheres, for uniform circular distributions of continuous simple 
poles in C and R^ immediate counter examples can be constructed. Third, in an attempt to derive the 
wanted uniqueness results, it is tempting to try to directly apply the standard theory of poles and residues 
associated with analytic function theory. For example, applying the residue theorem by taking a closed line 
integral around the unit disk immediately shows that X^s^i Mfe = however, additional progress quickly 
becomes difficult since, in order to make further progress, it is necessary to consider paths that extend into 
the interior region, where some form of analytic continuation inside the unit disk must be used, but this is, 
at best, questionable. These issues are addressed further in Appendix A. 

For simple poles the desired uniqueness theorem is: 

Theorem 2. If h{z) has the form specified by where is finite, < \zk\ < I and Zk' ^ Zk for k' ^ k, 
then h{z) = for all \z\ > I if and only if fi^ — for all k. 

Proof. Here /ifc = for all k trivially implies h{z) = for all z, so only the converse needs to be considered. 
Throughout the proof assume that < |zfe| < 1, z^:' 7^ for fc' 7^ fc and that \z\ > 1. The proof will be by 
contradiction, so assume to the contrary that a non-unique expansion exists where fik 7^ 0, but h{z) = for 
all |z| > 1. If /ife = occurs for any k in this expansion, then drop out these terms, reindex the /x^'s and 
reduce the value of Nk , so that ^ can be assumed to hold for all k without loss of generality. Using the 
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geometric series, each of the pole terras appearing in (jl4p can be reexpressed as 



- Zfe ^ 
n=0 



which has the same overall form as the power series 



71=1 



where a„ G 



(18) 



(19) 



As pointed out in many elementary treatments of complex variables, f{z) = for all z if and only if a„ = 
for all n > 1, which is a useful condition here. Substituting ^TE\i into allows h to be rewritten in the 
form: 



h{z) 



biZ , where b-t 



k=l 



(20) 



By assumption h{z) = holds for some given z^'s and /ifc's with Zk ^ z^i and /ife 7^ for fc = 1 2, 3, • • • Nk 
with A^fe > 0. Also as noted above, the factors bj occuring here are unique so bj = for j = 0, 1, 2, 3, • • • . 
From (j20p this condition on the bj 's can be rewritten in matrix form as 



G ^ = where G :— 



(1 


1 


1 • 






Z2 


■ 




4 


^2 


Z3 


■ 72 


4 


^2 


■ 

^3 


• z3 


V- 









(21) 



and where fi = (/xi, H2, /is, • • • , UNk)'^ ■ Consider the square matrix formed by the first Nk rows and Nk 
columns of G. This square matrix is the Vandermode matrix, which has a determinant that is well known 
to be nonzero when the Zk are distinct [3]. Thus the only solution to (|2ip is /i = 0, contrary to our original 
assumption, proving the uniqueness of the form given by (I14p . □ 

An analogous result can easily be shown for the logarithmic form given by (|10l) 



9{z) 



E 

fc=i 



Pk V'fe(^) 



where p G C is the source parameter and, as before, 

V'fc(z) = In 



{z - Zk) 



(22) 



(23) 



Here, since = implies ipk{z) = 0, so z 7^ and |z| > 1 > |zfc| > is assumed (as usual). Notice that from 
the discussion given in Appendix A, ipkiz) has no branch cuts for \z\ > 1, which means that it is analytic 
for |z| > 1 and thus that it has a proper power series representation: 



n—1 



(24) 
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Although (|24p can be obtained in various ways, it is easy to see that it is the correct series by simply 
comparing the derivative of the RHS of ([M]) with the series expansion of the derivative of the RHS of 
Substituting the series expansion into (|^ gives 



9i^) -EE 



fc=l n=l 



(25) 



and setting the various powers of z" to zero yields an equation set similar to ((21]) that must hold if g{z) = 
is to hold: 

/ Zl Z2 Z3 ■■■ ZATfc ^ 



G' p = where G' 



zf/2 zlll zlll ... z^^/2 
zll-i zllZ zll-i ... zlji 



\ 



(26) 



and where p = (pi, p2, Pa, 



As before, it is necessary to consider only the first TV^ rows of G', 



which can be reexpressed as a simple matrix product of the Vandermode matrix G and two other x TVfc 
matrices. If particular, if G' denotes this iV^ x A^fc matrix, then 



G' := N-^GX, where N 



/I 

2 

3 

Vo 








and X 



/zfe 

Z2 

Z3 

Vo 



\ 





(27) 



The condition G' p = thus becomes N^^GX p = 0, which immediately implies p = since all three matrices 
involved are invertible. This result can be restated as the desired uniqueness of logarithmic expansions given 
by the form (fTO|) : 



Theorem 3. If g{z) has the form specified by ilO\} . where Nk is finite, < \zk\ < 1 and Zk' ^ Zk for k' ^ k, 
then g{z) — for all \z\ > I if and only if p/j = for all k. 

Next consider the issue of the uniqueness of higher-order pole expansions of the form 



(28) 



where m is a positive integer, z^j, G C and the usual restrictions apply to z and z^. A series expansion for 
(z — Zfe)^™ = z^™(l — Zk/ z)^"^ can be found directly from the binomial series 



1 



(1 - 



m(m + 1) o m(m + l)(m + 2) 
1 + mx H -X^ H ^ rr^^ -X"^ 



2! 



3! 



mim + 1) . . . (m + n — 1) „ 
— 1 L 2 ia;" - 



(29) 



The general coefhcients here are related to the usual binomial coefficients. Introducing the special (and 
non-standard) symbol Sm,n for these signed coefficients and replacing x with Zk/z thus produces 



(Z - Zfe)™ z 



— E'^™"(~) where S'„ q := 1 and S„ 



n=0 



(rn + n — 1 ) ! 
n\{m — 1)! 



for n > . 



(30) 
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The expression of interest is thus 



1 _ c _ 

'^'-'i^)-i,j:^T.^>'^=^^ (31) 



z 

n=\ fc=l 

which leads to the foUowing matrix condition for the coefRcients: 

GXiy = (32) 

with V := (i^i, 1^2, 1^3, • ■ • , I'Nk)'^ ■ The solvabihty of this matrix equation directly yields a theorem expressing 
the uniqueness of higher-order pole representations: 

Theorem 4. For finite m > 0, if h^"^''{z) has the form specified by 128\) . where Nk is finite, < |zfc| < 1 
and Zk' 7^ Zk for k' ^ k, then = for all \z\ > I if and only if Vk = for all k. 

Here as m — > oo the theorem no longer holds since (z — Zfe)^'" — > for all \z\ > 1. Also Theorem 2] 
pertains only to poles of the same order. For very small Nk it is easy to show that expansions of mixed 
orders of poles are linearly independent by direct algebraic means. 

5 Uniqueness Results for M? 

Comparing ([151), dl) and (HT]) with shows that 



Re{h{z)} = Y. ^'^•^ (l^-^feT') 



fc=i 

where Dk = {—ak, —Pk)"^, which immediately shows that dipole expansions in are unique by Theorems[T] 
and[21 In a like fashion, it is clear that Theorem|3]in conjunction with Theorem[T](with a passive transforma- 
tion coefficient used as needed) also implies that higher order multipole expansions of order n are unique 
if they can be fully represented by two basis functions at each location. As noted in Section [21 in addition 
to the n — I case, this program can be carried out for n = 2 and 3, but for n > 3 there are more than two 
independent multipole basis functions for each Xk and so this correspondence cannot be one-to-one. 

First, it is useful to build on the multipole notation introduced in conjunction with ([TB|l and ([TT} by 
letting an i, i' , j or j' subscript preceeded by a comma denote d/dx when the subscript in question takes 
on the value 1 and to denote d/dy when it takes on the value 2 (in what follows it will be assumed that i, 
i' , j and j' always take on the values 1 to iV in and that similar partials are implied for R^, but only 
the case will be considered in this section). Multipole basis functions of any order can then be defined 
by expressions of the form 



II' ]...]' 1,1', 

where (as noted before) the number of subscripts corresponds to the order of the multipole. For n — 2 there 

ll^'i and Mi'Jj since obviously Ml^{ = M^^Jj 



are only two independent basis functions, say M!,^'^ and M!|^2 since obviously = M'i2 a-^d from Laplace's 



equation M2'°2 = — M!|^[. Likewise for rt = 3, if M!|*[ and M222 are selected as the two independent basis 
functions, then M!|^'^2 — ^^2221 M2'°2i — — M!|*!^ ^ and analogous relationships hold for permuted indices. 
For n = 4 it is easy to check that all of the other basis functions can be obtained from M.^^\ ^ , M2'°2 2 2 
and M^^\ ^ 2 • This shows the set of quadrupole basis functions for some given source point cannot be larger 
than the indicated set, but it does not show that the indicated sets are indeed independent. Showing the 
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actual independence of the indicated basis functions at a common source point follows from direct algebraic 
manipulation and involves first calculating them explicitly and then multiplying through by the common 
denominator, but this straightforward algebraic manipulation will not be done here. 

For n < 4, since there are only two independent basis multipole basis functions of order n at each source 
location, it is useful to introduce a more compact notation for these cases so that the resulting potential can 
be written more efficiently in terms of the above independent basis multipole functions 

i/(")(X)^^(4")Ai")+&i")Bi")) (33) 



where, for n = 1 
for n — 2 
and for n = 3 



k=l 



A^^^ = M^"' and B^^' = M^"' , (34) 
A^^^ = M}^{ and M^^'^ = M'"^ , (35) 
4^^=M''ii and Bf^ EE MtfJ^a ■ (36) 



Also in ([55]) the coefficients a^" and 6^" e M. Subsequently the superscript indicating the order is often 
omitted from a^,""* and fe^""*. 

Explicit expressions for aI"-* and bI^'' can be easily written: 



4^'-|:ln(|X-X.n^-^f^ , ^^-.^^ (37a) 



fel 



(37b) 



j^(2) ^ {X-Xkf - {y-Vkf g(2) ^ - Xk){v - Vk) 

^f^ ^^^'A i{y-ykf-{x-XkY] and Mf= ^^y'l']M x-xuf-{y~yk?]. (37c) 



\X-Xk 



The higher order complex poles corresponding to the expressions in p7p for n — 2 and n — i can easily be 
found, as in (jlSp . by multiplying the numerator and denominator of fik/iz — z^)" by {z* — z^)": 

p f f^k \ ^ ak[{x ~ Xk)^ - {y - yk)^] 2[3k{x ~ Xk){y ~ yk) .„„ . 

\{z-ZkYi \Z-Zk\^ \z~Zk\^ ^ ' 

f A^fc 1 ^k{x - Xk)[{x - Xkf - ijy - t/fc)^] /3fc(y - yfc)[3(a: - a:^)^ - {y - ykf] x 

Since \X — Xk\ ^ \z — Zk\, from ([37]) it thus follows that ([38|l can be immediately be rewritten as 

Re I , 2 1 - a,Af +/3fcBi') (39a) 

Clearly equations ([55]) imply an invertible passive coefficient mapping, so that Re {ft") (z)} ^ i?(")(X) 
holds. Thus Theorem 0] implies that multipole expansions of order three or less in are unique and these 
results can be summarized in the following formally as: 
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Theorem 5. In R^, for finite < n < 4, i/ H'^'^\X) has the form specified by i3S\) where Nk is finite, 
< \Xk\ < 1 and Xk' ^ Xk for k' ^ k, then = for all \X\ > 1 if and only if Ok ^ and bk = 

for all k. 

or less formally as 

Theorem 6. A finite expansion o/R^ points masses, point dipoles or point quadrupoles is unique. 

Since there are more than two multipole basis functions for each source location of an M'^ multipole of 
order four or higher, no simple correspondence exists between a given single higher order complex pole and 
the corresponding multipoles at the same location. Given that all of the uniqueness results considered so far 
rest on such a correspondence, it is clear that most of the readily obtainable results in M'^ have been found 
and multipoles of order n > 3 will thus not be considered. The next issue to be addressed is extending these 
R2 results to or even M^, for iV > 3. 

6 Uniqueness Results for 

This section addresses the question of point mass uniqueness by showing that uniqueness of point mass 
distributions in has direct consequences in . Towards that end the following simple lemma will prove 
to be useful: 

Lemma 1. In K^, for N > 1, an array of Nk distinct points forms at most 2 Nk x (Nk — 1) unique directions 
when all possible lines containing two or more array points are considered; moreover, it is possible to rotate 
coordinates so that all of the points are distinct when projected into the N — I dimensional hyperplane 
orthogonal to some preferred coordinate axis. 

Proof. Since two points determine a straight line, the maximin number of independent lines occurs when no 
three or more points are coUinear; i.e., chose any of the Nk points and a second distinct point. Factoring in 
the fact that any line determines two possible directions yields 2 Nk x {Nk — 1) for a maximum number of 
possible directions. (Besides the fact that many points may be collinear, many of the lines formed may be 
parallel, so the actual number might be much smaller than 2 Nk x {Nk — I)-) This shows the first part of 
the lemma. 

Next consider the second part of the lemma. Given a choice of coordinate origin, there is a continuous 
choice of preferred axis directions. Specifically, for a preferred coordinate axis (A^'th axis) choose a direction 
that does not coincide with any of the finite possible directions along which two or more points line up. Since 
no two points line up, their projection into the orthogonal hyperplane of the preferred direction is distinct 
and the lemma follows. □ 

This lemma leads to the following lemma: 
Lemma 2. // a point mass uniqueness counter-example exists in , then one exists in R^. 
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Proof. If a counter example exists in then from ^ it can be assumed, without loss of generality, that an 
Nk exists with nik 7^ for all Nk > k > 1 such that 

for all \X\ > 1 (where Xk' 7^ X^ for all k' 7^ fc). As before, without loss of generality, it is assumed that 
J2k=i = and Xk ^ for all k. Thus, for later convenience a set of point masses J2k=i = at the 
origin can be explicitly added to (|40|) . 

Let X = {x, y, z)'^ and Xk = {xk, Uk, Zk)'^ where there is no danger here of confusing z and Zk with 
the complex variables introduced earlier. Assume that the coordinates have been chosen in accord with 
Lemma [1] where the preferred direction has been taken to be along the z— axis so that in the x — y plane the 
Nk points are all distinct, so that {xk, yk, 0)"^ 7^ (xk' , yk', 0)^ for all k' 7^ k. Clearly a linear superpositions 
of potentials along the z— axis also obey the criteria of ((40l) so that 



nik ruk 

z=-L 

Observe that 



.=1 y\X-X,\ \X 



L L 

I dz .. f dz 

lim / — — = lim 2 



i^oo J y^{x - xk f + {y- VkY + {z~ ZkY J ^{x - XkY + {y- ykY + z'^ 

z=—L 2=0 



and that since 

L 



z=0 

it follows that 



[ , '^^ dz = In (l + ^Ja^ + L'^ \ - In a 
J \Ja} + z2 V ) 



L 



' ^21nr + V>ZE -21n 



z = -L 



+ z2 + z2 y I l + ^(6/L)2 + l 



where a = — x^)^ -\- {y ~ yk)'^ and 6 = -\/x^~T~y^. It thus foUows that 

/ 1 1 



Z — — 00 



\X-Xk\ \x\ 



dz ^2^k{X) (42) 



and thus that 

^ mk-^kiX) = with TOfe 7^ for aU fc. (43) 
fc=i 

□ 
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A comparison of (|43p and ([7]) allows one to one to conclude from Lemma [5] and Theorem IH] that 

Theorem 7. A finite expansion of R"^ points masses is unique, where the usual conditions are assumed to 
apply. 

Next consider the question of uniqueness for the dipole expansion given by ([5]): 

W{X) = J2 ^k-y {\X - . (44) 

fc=i 

It is possible to state and prove a dipole analog of Lemma [21 

Lemma 3. If a point dipole uniqueness counter-example exists in R'"^ then one exists in R^. 

Proof. Suppose, as in Lemma[2l to the contrary that an counterexample exists and thus that for some 
distinct Xk, that for some W given by the right hand side of W = Q for some \Dk\ ^ 0, for all k. From 
Lemma [T] there are only a finite number of directions parallel to the lines determined by two or more of the 
Xfe's. To this finite collection of directions add the vectors {-Dfe}^-^ and then select a preferred z coordinate 
direction that is different from all of these directions and denote this direction by k. Since \k ■ Dk\ < \Dk\ for 
all k and, by construction, not only does {xk, yk, 0)'^ 7^ {xk', yk', 0)"^ for all k' ^ k hold, as in the proof of 
Lemma [21 but the projection of Dk on the x — y plane is nonzero. Let Z?^ denote this projection. Further, 
as before, the case J2k=i -D/c 7^ presents no real difficulties and {xk, yk, 0)^ 7^ can also be assumed. Then 
taking the integral along the z— axis as in (|42)) gives [after adding an analogous term at the origin to the one 
added to dH])] 

dz =2-^^kiX) (45a) 



Z— — 00 






1 


1 


\x 


-Xk\ 


\x\ 




1 


1 


\x 


-Xk\ 


\x\ 




1 


1 




-Xk\ 


\x\ 



X 



dz =2|-*fe(X) (45b) 
dy 

dz =0 (45c) 



Using (|45p in (|44ll allows the R^ counterexample to be restated as 

/ WiX)dz = 2^ ^f^-V*fc(x, y) = (46) 

where ^| ^ for aU k. □ 

Since Theorem [6l shows that a counterexample of the form specified by (|46l) cannot exist, the following 
theorem is an immediate consequence of Lemma [31 
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Theorem 8. A finite expansion o/M^ point dipoles is unique. 

Clearly, one would expect the results of Lemma [5] and thus Theorems [7] and [5] to generalize to for 
> 3, with only the technical difficulty of explicitly obtaining the N dimensional integral analogs of 
and (|45p to stand in the way; however, two other remaining issues of far greater practical importance are 
less clear: (1) Proving quadrupole and other multipolc uniqueness results in . (2) Proving expansions 
of mixed types of point sources are unique. It is unclear how the first issue should be approached, but an 
approach to the second issue can be based on the theorems of Section [3l As before, results in the complex 
plane will be the starting point. 



7 Uniqueness for Point Sources of Mixed Type 

This section deals with uniqueness results for expansions of mixed type. Only the complex setting will be 
considered here even though results for mixed types in C can be directly extended to and mixed type 
results. There is a pragmatic reason for considering only the complex setting: Uniqueness results for mixed 
types of expansions are very difficult to prove and, in the end, given that there is only a limited amount of 
success here in studying the complex case, consideration of the real case is irrelevant. 

Thus two particular kinds of mixed type analytic expansions will be of primary interest here: (1) Expan- 
sions consisting of both simple pole terms and second order pole terms. (2) Expansions consisting of simple 
poles and logarithmic basis functions. 

First, observe that as far as uniqueness results for expansions of mixed type are concerned, it only 
necessary to consider the general case where there are (possibly) a like number of either simple poles and 
logarithmic basis functions that are located at the same point or, in a like manner, to consider only the case 
consisting of simple poles and second order poles that are located at the same point. Thus, without loss 
of generality, for the simple pole and second order pole case consider the following expression for poles of 
mixed type: 

2) ^ V + E (47) 



1 . ^ r. 



For m = 2, jSQl) gives 

7 ^ = ^> '(" + !)—, (48) 

which can be combined with (|18|) to yield 



The RHS of (SQl) can be rewritten as 



z '■ — ' '■ — ' z" z^ '■ — ' '■ — ' Z'' 

k=\ n=0 k=l n=0 



Z'^ M-^ M-^ nZi 



z-h^-'Hz)=j:,,+j:^^Eii +E-^-E^- (50) 

k=l k=l n=\ k=l n=l 

Setting the successive powers of z^" individually to zero on the RHS of (j50|l (and simply ignoring the 
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condition X^fe^i /^fe = ^) gives the following set of equations that must hold when z-/i(^'^^(z) = 0: 





Z2 


Z3 ■ 




1 


1 


1 


1 \ 


4 


^2 


Z3 




2zi 


2Z2 


2Z3 


2zAr, 


^3 
^1 


^3 
Z2 


Z3 














^2 


^3 




NkZ^"- 


NkZ2 


N1.-1 
NkZ^'^ 


JVfc -1 




"fe + i 
^2 


^3 




{Nk + l)z'^' 


(Nk + l)z^'' 


(iVfc + l)z^'= • 





/ ^1 \ 

M2 
1^3 



= 



(51) 



The first 2Nk rows of this matrix equation set can immediately be rewritten in the following block matrix 
form: 



GX 


NG 




(N + 7VfeI)GX"" 



0. 



(52) 



Here I is the Nk x Nk identity matrix and the other matrices in ([5^ have been previously defined. (|52p is 
equivalent to the two coupled equations sets 



GX /i + NGz/ 







GX^'=+V+ (N + iVfcI)GX"V = 
Here the first of these equations uniquely determines fi in terms of v. 

= -X-^G-^NGi^ . 



(53) 
(54) 

(55) 



When this is substituted into the second matrix equation and the result is rearranged slightly the following 
matrix equation for results 



[GX^'=(iVfcI- G-^NG) + NGX"']j/ = , 
which can be multiplied from the left by X^^'-'G^^ to give 

[(ArfcI-G-iNG) + X""'"G-iNGX"'=]i/ = 0. 



(56) 



(57) 



Here ([55)1 immediately implies that if the matrix [(TV^I — G ^NG) + X ''G ^NGX is invertible then 
V = Q must hold. In turn, this matrix is invertible if the following matrix is invertible: 



which can be rewritten as 



C := [TVfcl-N + GX-^^G-iNGX^'G-i], 
= iVfcl - N + U-^NU, where U := GX"*" G'^ 



(58) 
(59) 



Showing that C is invertible is not as easy at it might first appear. Thus, for example, first observe that C 
can be rewritten as C = A + B, where A := Nkl — N and B :— U~^NU. Obviously A is positive definite; 
furthermore, it is clear that the sum of two positive definite matrices is positive definite, so if it can be shown 
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that B is positive definite, then C will be positive definite and thus invertible. Next observe that it is trivial 
to find the eigenvectors of B and that the corresponding eigenvalues are given by the diagonal elements of 
N. As one can quickly convince him or herself, this, however, does not imply that B is positive definite, 
since, among other things, the eigenvectors are not orthogonal to each other (clearly B is not normal since 
BB^ 7^ B^B and normality is assumed in relevant theorems of interest). Given that the author's attempts 
at proving that C is invertible have not met with success, the issue is open. 

Here the question naturally arises as to what the analog of the C matrix would have been if the set of n 
rows starting at row mn + 1 had been taken instead of at row n + 1. In this case the system of equations 
that result are 

Cmiy = , (60) 

where: 

C„, := miVfel - G-^NG-i + x"""' G-iNGX""' . (61) 

Here it is also natural to also raise the question about what the analog of the C is when logarithmic basis 
functions are included. Thus consider the conditions that must hold if (p{z) = h''^'^'^\z) = 0, where 



fe=l 



Since 



dipk{z) 



Hk 



{z - Zk) 



dz 



Z- Zk 



■E- 



+1 



{z - ZkY 

oo 

-E 



oo ,n-l 



n=2 



it is obvious that the correct power series expansion for 'i()k{z) is 

Mz) = Y.-ii 

n z" 

ri=l 

Substituting this expansion along with 



(62) 



(63) 



(64) 



Zk 



n=0 



n=l 



(65a) 



and 



yields 



n=2 



(n-l)4" 



{z - ZkY 

/^^^^^^^H^) = E[E?^I + E^ 



oc n-1 



+ E 



{n-l)zl 



-fk 



fc=l n=l n=l 

Breaking out the n = 1 terms separately and reindexing yields: 



n=2 



h^^'-'^\z) = W^PkZk+l^k) + W-nT. 



k=l 



Z — ' Z' 

n=l fe=l 



PkZk 



n+1 



(n + 1) 



+ z^lj,k+nz]^ ^Vk 



(65b) 



(66) 



(67) 
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Introducing p'^, = ZkPk and v'^, = [ikjzk yields 



n=l fc=l '-^ ' 



(68) 



Ignoring the first term on the RHS of (|68p and setting the other powers of 1/z to zero yields an equation set 
whose first 3n rows can be written in block matrix form as: 



(N- 




iG 


GX 


NG 


[N + (TVfe ^ 


-1)] 


-il)GX"'^ 


GX"*" 


(N 


-f iVa)GX"'= 


^ [N + (2iVfe - 


HI)] 


-il)GX"^'' 


GX^"'' 


(N4 


2iVfeI)GX'"'' 



= 0. 



(69) 



When there are no second order pole terms, (j69p can be rewritten as 

P' \ 



G 


(N + I)GX 


GX"^ 


[N+(iVfe + l)I]GX"'' 



= 0, 



which, by following analogous steps used to obtain ()57p . can be rewritten as 

[(iVa-G-iNG)+X""'"G-iNGX"']Ai = 



(70) 



(71) 



Hence combined logarithmic and simple pole basis functions are independent if the C matrix introduced 
earlier is nonsingular. 
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Appendix A 



Branch Cut and Analytic Continuation Side Issues in C 



This appendix addresses two distinct questions (or misperceptions) that some readers may have: (1) Since 
the logarithmic point source basis functions tpk [as given by ^] are used, and logarithmic functions generally 
have brach cuts in C, does ipk have problematic branch-cuts and, if not, why not? (2) When only simple poles 
may be present in the interior region, since it seems natural to assume that the analytic continuation of a zero 
function is always a zero function, can one extend analytic continuation into the interior of the unit disk along 
various paths while assuming that f[z) = still holds until there are areas of overlap, so that one can simply 
integrate around each pole separately and then directly prove the desired complex plane uniqueness results 
for simple poles from a direct application of the residue theorem? This analytic continuation procedure may 
be tempting, because for a finite collection of potential simple poles in the interior. Theorem 2 guarantees 
that if f{z) = holds in the exterior region, then = 0, so this analytic continuation procedure seems to 
work for this case. As far as counter examples go, one need only consider a simple circle with a uniform 
density constant simple pole strength and a compensating interior pole, but what about a denumerable set 
of separated simple poles, which is the case of prime interest here? The reader who is not bothered by this 
last sort of question may simply skip the second part of this appendix that deals with this issue. 



Branch-cut Issues 

First consider branch-cut related issues in the exterior of a unit disk, when all the source points reside 
inside the unit disk. Although simple and higher order poles obviously do not have such branch cuts, at first 
glance it might appear that 4'k{z) given by ^ does have branch cuts: 

Mz) ■■= In — In - (A-1) 

Z~ Zk z 

since, for example, 

In z ^\nr + i6 . (A-2) 

The truth of the matter, however, is not so straightforward and it turns out that there are no branch cuts 
in the exterior of the unit complex disk for ipkiz). To see this geometrically, first let £k ■— \z — Zk\ and 
let (j)k equal to the angle between the positive x-axis direction and the vector parallel to the line segment 
connecting Zk to z [which is to say, the angle between the line segment connecting the points {xk, Uk)'^ 
and (x, Uk)'^ and the line segment connecting [xk, yk)"^ and {x, y)'^]- Then 

Mz) := In f (A-3) 

and when one draws a plane figure displaying the various relevant lines and angles, it is perfectly obvious 
that, for any choice of z^^, and (j)k are equal in value at two places as z traces out a closed loop around 
the unit disk (keeping in mind that \z\ > 1). In fact, from the geometry of this figure it is clean that 
tt/2 > \9 ^ (j)k\ must always hold and, hence ipkiz), for each k. is uniquely defined for z\ > 1. 
This branch-cut issue can also be settled using analysis by considering the form 

Mz) = In := In \ . (A-4) 

[Z - Zk) (l~Zk/Z) 
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and noting that since Re{l — Zk/z} > (which follows immediately from 1 > \zk/z\ > |Re{zfe/z}| ), the 
absolute value of the argument of 1/(1 — Zk/z) is less than n/2 from elementary properties of the arctan 
function. 

Notice that while there are no branch-cuts for 4'k for \z\ > 1, for |z| < 1 it is clear that there is a branch 
cut connecting Zk and the origin. An expansion of the form (jlOp thus has a collection of branch cuts that 
form a star-like pattern in the interior of the unit disk. 

Finally, since a linear superposition of analytic function is analytic the use of {ipk}k=i a set of basis 
functions, as in (|10p entails no branch-cut interpretational issues at all. It is clear, however, that one cannot 
decide to treat the RHS of a t/jk fit, or expansion such as (|10p . as a single logarithmic function using the 
standard properties for combining logarithms-when a linear combination of V'fc's are combined into a single 
composite logarithmic function brach-cut issues are arbitrarily introduced. (Since the argument of the log 
of some term is assumed to be between and 2tt, even absorbing the constant /ifc into ipk by itself causes 
problems and introduces unwanted restrictions, because the complex part of iik'tpk is not similarly restricted.) 

Analytic Continuation Issues 

Next, consider the analytic continuation issues, which are perhaps best addressed by consideration of 
counter-examples. It may seem reasonable to argue that continuation of the zero function is a special case, 
especially when it done on either side of a neighborhood where it is known that, at most, one simple pole 
resides and it would seem that what happens outside the greater region under consideration does not matter. 
Thus, suppose that an expansion of the form is being considered, and that the poles are ordered such 
that zi is the location of the simple pole that is closest to the analytic region (|z| > 1). Consider a "fit" to 
the function f{z) — 0. Then, if one analytically continues the zero function into the unit disk close to zi 
along a path to one side of this pole, it would seem obvious that f{z) = over this entire region. On the 
other hand, if one starts from the same region and does the same thing on the other side of the pole then 
in the region of overlap between these two analytic continuations, it is obvious that f{z) ~ 0. From this, 
one can conclude that the simple pole necessarily has a weight of zero: [lu = 0. Proceeding in a like fashion 
to each succeeding pole, one could thus argue that /(z) = for \z\ > 1 implies that fik — for all fc, which 
is the desired result. Although, at each step, since there are two analytic regions of continuation that are 
not simply connected, and it is known that analytic continuation into regions that are not simply connected 
are problematic; the simple pole in question produces no branch cuts in either of these two regions and, 
furthermore, the f(z) always agrees in the region of overlap. No function is simpler than the zero function 
and one might argue that in this special case of analytic continuation there is no problem, at least in this 
particular instance, and, when used in this way, it leads to the right result. 

To see what is wrong with this argument, consider the following. It will only be necessary to show 
that one counter example exists (here some liberty will be taken and it will merely be shown that some 
counter-example is very likely to exist, without explicitly constructing the explicit counter-example itself). 
Let {sj}fLi be a sequence of points in the interior of the unit circle (excluding the origin). Consider the set 
of corresponding points in the exterior of the unit circle given by tfc = 1/s* . Then it is well known that the 
sequence of values {f{tj)}j^^ completely characterizes the analytic function f{z). Consider a DIDACKS fit 
of the form 

^(z)=f;^^, (A-5) 

which always exist for N < oo since the associated linear system was always shown to be solvable in [7]. It 
does not seem unreasonable to assume that for some sequence of points {t^} and some choice of analytic 
function f{z) that this system remains solvable as iV ^ oo. [For example, it is clearly perfectly acceptable 
to choose an f{z) that has the form (jA-5[) itself, with some appropriate choice of {aj}°^^ and {sj}°^i, so 
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long as \aj\ ^ (in the argument that foUows, it is important that the sum in (|A-5p contains an infinite 
number of simple poles).] 

Now, f{z) for \z\ > 1 could have been represented by its values at some other infinite sequence of points, 
say {/(Pfc)}fe°=i; where pk ^ tj for all j and k. Here let Zk = and assume that a fit of the form 
with n — oo also exists. Then, for \z\ > 1: 

oo oo 

Since by construction, aj ^ for all j and 7^ for all fc, it is clear that the strategy of analytically 
continuing the zero function so as to encompass an isolated simple pole must fail if any representation of 
f{z) = exists that has the form (jA-6|) : moreover, for this to occur it is only necessary that for some aj 7^ 
for all j and some fik ^ for all fc, that sj and pk exist such that 

oc oc 

E^ = E^' (A-7) 

^ — ' Z — Sk ^ — ' Z — Zk 

k=l ^ k=l " 

where (for all j and fc) < \sk\ < 1, < \pk\ < 1 and Sj ^ pk- It seems most probable that two such 
bounded sequences of simple poles exist. 
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